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ABSTRACT 

This  paper  considers  the  problem  of  stabilizing  a  control  system  governed 
by  a  combination  of  partial  and  ordinary  differential  equations.  The  partial 
differential  equations  govern  the  evolution  of  the  system  in  the  interior  of 
some  spatial  domain,  the  ordinary  differential  equations  describe  the  evolu¬ 
tion  of  the  boundary  data;  the  control  enters  through  the  boundary  ordinary 
differential  equations  in  a  bilinear  fashion.  We  provide  sufficient  condi¬ 
tions  for  feedback  stabilization  of  such  ^hybrid*  systems.  Two  examples  to 
wave  equations  with  dynamic  boundary  conditions  are  provided. 
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SIGNIFICANCE  AND  EXPLANATION 


In  the  control  of  distributed  parameter  systems  i.e.,  systems  governed  by 
partial  differential  equations,  controls  often  times  can  be  applied  only  on 
the  boundary  of  some  spatial  domain.  In  this  paper  we  consider  the  case  when 
the  controls  enter  the  boundary  conditions  as  bilinear  (in  the  control  and 
state)  ordinary  differential  equations.  We  show  how  we  can  synthesize  feed¬ 
back  controls  that  will  stabilize  systems  whose  uncontrolled  motion  is 
critically  stable. 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MRC,  and  not  with  the  authors  of  this  report. 


FEEDBACK  STABILIZATION  OF  "HYBRID"  BILINEAR  SYSTEMS 


♦  * 

M.  Slemrod  and  E.  L.  Rogers 

0.  Introduction 

In  a  recent  paper  [1]  Burns  and  Cliff  formulated  a  model  of  a  "hybrid"  system, 
i.e.,  a  mixed  system  of  partial  and  ordinary  differential  equations,  in  which  the 
control  enters  only  within  the  context  of  the  ordinary  differential  equations. 
Interest  in  such  systems  can  be  motivated  by  problems  in  structures  in  which  the 
control  dynamics  take  place  only  in  the  boundary  conditions  of  a  distributed 
parameter  system.  In  this  paper  we  consider  two  problems  (motivated  by  the  example 
of  Burns  and  Cliff)  where  the  control  enters  the  boundary  conditions  in  a  bilinear 
fashion.  Our  goal  is  to  synthesize  feedback  controls  which  will  stabilize  the 
originally  critically  stable  systems,  i.e.,  we  wish  to  find  a  feedback  controls  so 
that  the  states  of  the  feedback  systems  approach  the  zero  state  as  t  *  a*. 

The  main  tool  of  our  analysis  will  be  "hyperbolic"  stabilization  theory  of  Ball 
and  Slemrod  [2],  [3].  In  fact  the  theory  of  [2],  [3]  may  be  readily  applied  to  our 
problem. 

The  paper  is  divided  into  five  parts.  Section  1  recalls  the  results  on 
feedback  stabilization  of  [2],  [3].  Section  2  introduces  the  first  hybrid  system. 
Section  3  shows  the  feedback  stabilizability  of  the  first  hybrid  system.  Similarly 
Section  4  discusses  the  second  hybrid  system  and  Section  5  proves  feedback 
stabilizability  for  it. 
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1.  Feedback  Stabilization 


Consider  the  initial  value  problem 

u(t)  *  Au(t)  +  p(t)  u( t) /  (P) 

u(0)  -  uQ  e  X, 

where  A  is  the  infinitesimal  generator  of  a  Cq  semigroup  of  contractions 

e^  on  a  real  Hilbert  space  X  ,  8  is  a  bounded  linear  operator  X  ♦  X, 

and  p(t)  is  a  real  valued  control.  X  is  endowed  with  inner  product 

<*/•>  • 
x 

Definition  System  (P)  is  stabilized  (weakly  stabilizable)  if  there 
exists  a  continuous  feedback  control  p:X  ♦  R  such  that  ( P)  with  p(t)  * 
p(u(t))  satisfies  the  following  properties. 

(i)  For  each  uQ  there  exists  a  unique  weak  solution  of  u(t;ug) 
defined  for  all  ten,  of  (P). 

(ii)  {0}  is  a  stable  equilibrium  of  (P). 

(iii)  u(t,Ug)  ♦  0  (u(t,Ug)  ♦  0  weakly)  as  t  ♦  in  X  for  all 

Ug  e  X. 

The  natural  approach  to  the  stabilization  problem  is  to  differentiate 
2 

lu(t)lx  ■  <  u(t),u(t)  >x  along  trajectories  of  (P).  In  this  manner  we 
obtain 

^-lu(t)l^  -  2  <  Au(t),u(t)>  +  2p(t)<u(t),8u(t)>  , 

at  least  formally.  An  obvious  choice  of  feedback  control  ( though  not  the  only 
one)  is 

p(u)  ”  -<u,8u>x, 

since  this  control  yields  the  "dissipative  energy  inequality" 

~  lu(t)*x  <  -2<u(t),8u(t)>x  . 
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(Note  that  e 


a  contraction  means  <Au(t) ,u(t)>^  <  0,  for  u(t)  e  D(A))« 
So  formally  this  choice  of  control  p(u)  yields  a  feedback  system  of  that 
form 

u(t)  =  Au(t)  -  <u( t ) , 8u( t ) >x  u(t).  (F) 

For  the  purposes  of  this  paper  the  following  theorem  of  Ball  and  Slemrod 
[3]  will  be  needed. 


Theorem  1.  If  8 :  X  ♦  X  is  compact  and 

<eAS,  8eAtiJ»x  =  0  for  all  t  e  R+  ==>  ^  =  0  (C) 

then  (P )  is  weakly  stabilizable. 

In  some  aplications  it  will  be  convenient  to  work  with  second  order 
"hyperbolic"  systems  of  the  form 


y(t)  +  Ay(t)  +  p(t)By(t)  =  0,  (P‘) 

y(0)  =  yQ  e  Hftf  y(0)  =  e  H. 

Here  H  is  a  real  Hilbert  space  with  inner  product  <*,•>  and  norm 

H 

2 

1*1  =  <•,•>.  A  is  a  densely  defined  positive  self-adjoint  operator  on 

H  H 

H  such  that  A  *  is  everywhere  defined  and  compact.  We  suppose  the  eigen- 
2 

values  A  of  A,  n  =  1,2,***,  0  <  A.  <  \  <  ••• ,  are  simple.  We  denote 

n 

the  corresponding  sequence  of  eigenfunctions  by  },  Let  H.  =  D(A  ^  ). 

n  ** 

Ha  form  a  Hilbert  space  under  the  inner  product 

*  V-,  * 

<y.y  >A  *  <A  '2y,A  '2  y  > . 

Denote  I • I  as  the  norm  of  H. .  We  assume  B  is  a  bounded  linear  map  :  H 

A  A  A 

Again  the  function  p(t)  is  a  scalar  real  valued  control. 

If  we  set 


y(  t) 

u(t)  =  (  .  ), 

y(t) 


I 

0 


). 


8  - 


D(A)  =  D( A )  *  H  , 
A 


♦  H. 
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x  -  *  h  with  <(y,y)/(y  ,y  )>x  =  <y,y  >H  +  <y,y  >  . 

A 

we  see  A, 8  are  required  for  problem  ( P)  above.  We  can  now  state  t'-.e  theorem 
of  Ball  and  Slemrod  given  above  in  a  second  order  context  (see  [2)  for  a 
generalization ) . 


Theorem  2.  Suppose  A  and  B  are  as  above.  Assume  B  :  +  H  is  compact 

and 

(HI)  <»♦„,♦,>„  *  0  for  k  =  1 , 2 ,  •  •  • 

X  X  H 

(H2)  A  ±  X  *  2X  unless  m  =  n  =  k  and  the  +  sign  is  taken, 

m  n  k  si 

or  both  <B$  ,$  >„  and  <B$  ,  ♦  >  are  zero, 
in  n  H  n  m  H 


Then  the  feedback  system 


with 


y(t)  +  Ay(t)  +  <By(t) ,y(t) )  By(t)  *  0 

H 


(F*> 


p(y(t),y(t))  *  <By(t) ,y(t)> 

H 

and  the  initial  data  y(0)  =  yQ  e  HA,  y(0)  =  y^  e  H  possesses  a  unique 
globally  defined  weak  solution  (y,y)  e  C((0,»)?HA  x  H)  and  (y(t),y(t))  + 
(0,0)  weakly  in  X  =  Hft  x  h  as  t  ♦  08  and  (P* )  is  weakly  stabilizable. 


Proof.  B  :  ♦  H  compact  implies  8  :  X  ♦  X  is  compact  and  (HI),  (H2) 

imply  (C).  See  [2]  for  details. 
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2.  Hybrid  system  1. 

Motivated  by  the  example  in  [1]  we  consider  the  string-mass  system  shown 
in  Figure  1.  The  string  has  length  1,  constant  linear  density  o,  and  is 
under  constant  tension  T.  The  purpose  of  the  device  at  the  right  end  is  to 
maintain  the  tension,  however  it  is  idealized  so  as  to  provide  no  impedance  to 
the  vertical  motion  at  the  end. 

The  vertical  motion  of  the  string  is  assumed  to  satisfy  the  linear  wave 
equation 

ztt(t,x)  =  «\x(t,x).  0  <  x  <  1,  (2.1) 

2 

where  a  =  T/a  >  0.  The  motion  of  the  right  hand  end  of  the  string  is 
governed  by  the  balance  of  linear  momentum.  The  relevant  forces  here  are 
acceleration  of  the  point  mass  m,  the  tensile  form  on  the  string  Tz  (t,1), 
and  external  forces.  We  assume  we  can  impose  an  external  force  at  the  right 
end  in  the  bilinear  form  p(t)z(t,1).  In  the  absence  of  other  external  forces 
the  equation  of  balance  of  linear  momentum  of  the  mass  m  is 

mz  (t,1)  =  -p(t)z(t,1)  -  Tzx(t,1).  (2.2) 

The  left  hand  is  assumed  fixed  at  x  =  0  so  that 

z(t,0)  =  0.  (2.3) 

The  initial  conditions  are 

z(0,x)  =  f(x),  zt(0,x)  =  g(x) ,  0  <  x  <  1.  (2.4) 
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3.  Abstract  formulation  of  the  hybrid  system  1  and  feedback  stabilization. 
For  the  hybrid  system  (2.1) -(2. 4)  set 


z(t,x) 

y(t)  =  (  }  where  z  (t)  =  z(t,1). 

z  1  ( t ) 


Also  define  the  operators 


f  -«2  i- 


1  £l_ 

m  dx 


0  0 


1  ^  ' 


where  D ( A )  =  {{z,z^)  e  H  (0,1)  x  r  :  z(o)  =  0,  z(1)  =  z  }.  Set 

2 

H  =  L  (0,1)  x  R  where  H  is  a  Hilbert  space  endowed  with  the  inner  product 

<(v,w),  (v  , w  )>  =  “  J1  v( x ) v  (x)dx  +  —  w  w  . 

H  2  '  t 

a  0 

Clearly  D ( A)  is  dense  in  H.  Of  course  we  see  that  now  (2.1)-(2.4)  has  the 
form  (1.1). 


Lemma  1 .  A  is  a  positive  definite,  self-adjoint  operator  on  H,  with  bounded 
inverse  A-1, 


*  * 

Proof .  We  first  prove  A  is  symmetric,  i.e.,  <A(v,w),(v  ,w  )>  = 

H 

*  ★  *  *  A 

<(v,w),A(v  ,w  )>  for  all  (v,w),(v  ,w  )  e  D ( A )  where  D( A)  c  d(A  ).  To 

fi 


see  this  simply  compute 


<A( v,w) , ( v  ,w  )>  =  Z1  v'(x)v*  (x)dx 

H  0 

*  * 

=  <(v,w),A(v  , w  )>„ 


(3.1  ) 


Furthermore  since  for  (v,w)  e  D( A )  we  have  v(x)  =  Jx  v'(x)dx  <  (/1v'(x)2dx) 


and  hence  the  inequality 


sup  |v(x)|2  +  /1v(x)2dx  <  const.  J1  v*(x)^dx. 
0<x<1  0  0 


(3.2) 
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Thus  from  (3.1)  we  have 


2 

<A(v,w) , (v,w)>  >  const. I (v,w)l  ,  const.  >  0,  (3.3) 

H  n 

i.e.,  the  operator  A  is  accretive.  By  the  Schwarz  inequality  jo  have 

2 

lA(v,w)l  I  (v,w)  I  >  const.  Mv,w)l 
Ha  H 

and  hence 


*A(v,w)l  >  const.  I(v,w)l  .  (3.4) 

H  H 

Thus  A  possesses  a  bounded  inverse  on  H.  Hence  R(A)  =  H  and  a  well  known 
result  (e.g.  see  Yosida  [4;  p.  199])  shows  D(A*)  c  D( A }  and  A  is  self  adjoint. 
Inequality  (3.3)  shows  A  is  positive  definite  □ 


Lemma  2.  A 


H  +  H  is  compact. 


Proof.  Let  {(v  ,w  )}  be  a  bounded  seauence  in  H  and  set  A  *(v  .w  )  = 
n  n  n  n 

(z  ,z.  ),  n  =  1,2, •••  .  Since  A-1  is  bounded  l(z  ,z.  )l„  <  const,  for 

14  111  n  in  H 

all  n.  Also  since 


0  ,2  T  dz  (1) 

*2  d  i  t  /  x  ^  n 

-a  — ~  z  (x)  =  v  (x),  -  - - 

.  2  n  n  m  dx 

dx 

Multiplication  of  the  first  equation  in  (3.5)  by  z 


yields 


=  w  .  (3.5) 

n 

and  integration  by  parts 


®2  /1z'(x)2dx  =  a2z  (1)z'(1)  +  J1  z  (x)v  (x)dx 
0  n  n  g  n  n 

»  a2  S  z  (1)w  +  f1  z  (x)v  (x)dx 

*  n  n  g  n  n 

*  a2  <(z  ,z  ),(v  ,w  )> 

n  In  n  n 

<  a2l(z  ,z  )l  l(v  ,w  )l  <  const, 
n  In  H  n  n  H 
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Since  z'^(x)dx  >  /2  z^(x)dx  by  (3.2)  we  know  {z  }  lies  in  a  bounded  set 
0  n  0 

of  (0, 1 )  •  Since  the  injection  of  ^(0,1)  into  L2(0,1)  is  compact  {z^} 
possesses  a  convergent  subsequence  in  L^(0,1).  Since  {z1n}  a  bounded 
sequence  in  R  it  certainly  contains  a  bounded  subsequence.  Hence  A  *  maps 
bounded  sets  of  H  into  precompact  sets  of  H  and  hence  is  compact.  0 


2  2 

Lemma  3.  The  eigenvalues  of  A  are  given  by  a  sequence 

where  {X  }  are  the  increasing  positive  roots  of 
n 

— *  tan(^)  t  {(sinX  x,  sinX  )} 
muA  o  n  n 

are  the  associated  eigenvectors. 


Proof.  If  is  an  eigenvalue  of  A  ,  then 


-a 


2  d^z 

dx2 


,2 

A  z 


■  »<«<■-  i^r1 


X2z(1), 


and  z(0)  ■  0,  where  (z(x),z(1))  is  the  associated  eigenvector.  Thus  z(x) 

rX> 


sinf— )x  and 
ct 


x  r  Xt 


□ 


Lemma  4.  B  is  a  bounded  linear  self-adjoint  operator  :  H  ♦  H. 


Proof.  I B( v ,w ) I ..  =  1(0,— )l„  <  const.  I(v,w)l„  so  B  is  bounded.  Also 

1  H  m  H  H 

*  *  i  *  *  * 

<B( v |W) » ( v  ,w  )>  =  —  ww  =  <(v,w),B(v  ,w  )>  so  B  is  symmetric  and 

H  i  n 

bounded,  hence  self-adjoint.  □ 


Lemma  5.  <B(sinX  x,  sinX  ),  (sinX  x,  sinX  )  >..  t  o  for  n  =  1,2.' 
n  n  n  n  h 


. 

Proof.  <B(sinX  x,  sinX  ),(sinX  x,  sinX  ) )..  =  —  X  *  0. 
-  n  n  n  n  H  x  n 
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r 


Lemma  6.  For  the  eigenvalues  of  A  we  know  1  ±1  *  2X,  unless  m=n=k 

"  in  n  K 

and  the  +  sign  is  taken. 


Proof.  Assume  X  ±  X  =  2X.  .  Divide  by  a  and  take  tan  of  both  sides. 
-  m  n  k 


We  then  have 


X  X  2X 

tan(-~  ±  -£)  -  tan(— ^ ) 


a  a  '  'a 

which  in  turn  implies  by  the  usual  trigonometric  identities 

XX  X. 


2  tan  — 
a 


tan(-“)  ±  tan(~) 

1  * tan(;r)  tan(^r)  1  - tan2(^r) 


We  now  use  the  definition  of  X  (Lemma  3)  to  assert 

n 


or 


s  <C  *  2(^  C 

1  *  (Ir)2^"1x”1 

ma  m  n 


1  -  (— )2x"2 

ma  k 


(Xn  *  V 

W  <=>* 


2X. 


X2  *  >2 

k  ma 


(3.7) 


We  now  have  two  cases  to  consider. 


Case  1.  X  +  X  ■  2X.  . 

-  m  n  k 

2 

In  this  case  (3.7)  implies  XX  *  X  .  But  since 

n  n  K 

2 

see  (X  -  X  )  =0  and  X  »  X  =*  X,  . 

m  n  m  n  k 


X2  +  X2  +  2X  X 
m  n  m  n 


Case  2. 
In  this 

or 


X  -  X 
m  n 

X  -  X  * 
n  m 

X  X  +  ( 
m  n 


2 


and  (3.4) 


k  ma 


implies 

2 
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2 

XX  =  -A,  which  contradicts  the  positivity  of  {A  }. 
m  n  k  n 

This  completes  the  proof.  □ 

Theorem  3,  The  feedback  system  (2.1),  (2.2)  with 

z  <t, 1 )z(t,1 ) 
p(t)  „  -£ - - - 

possesses  a  unique  globally  defined  weak  solution  for  initial  data  (f,f(1))  @ 
Hft,  (g,g(1))  e  H,  ((y,y-,)»  (YiY^))  +  ( (0,0 ),  (0,0) )  weakly  in  *  H  as 
t  ♦  °°,  and  (2.1),  (2.2)  is  weakly  stabilizable  in  *  H. 

Proof.  Lemmas  1-6  show  the  hypotheses  of  Theorem  2  are  satisfied.  □ 


4.  The  hybrid  system  2. 

“J 

Suppose  G  is  a  domain  in  R  filled  with  a  compressible  fluid  which  is 
at  rest  except  for  accoustic  wave  motion.  If  $(x,y,z,t)  is  *■  n  velocity 
potential,  so  that  is  the  particle  velocity,  then  linearized  theory  says 

that  $  satisfies  the  wave  equation 

=  c2A$  in  G,  (4.1) 

where  c  is  the  speed  of  sound  in  the  medium,  is  the  pressure  distribu¬ 

tion  (  *  w(x,y,z,t)).  Now  suppose  that  the  boundary  of  G  3G  possesses  a 
non-rigid  section  which  is  subject  to  small  oscillations.  We  assume  that 

each  point  on  r  reacts  to  the  excess  pressure  of  the  accoustic  wave  like  a 
harmonic  oscillator.  We  assume  also  that  different  parts  of  the  boundary  do 
not  influence  each  other,  so  that  r  is  locally  reacting.  Then  the  normal 
displacement  6  of  the  boundary  I*  into  G  satisfies  an  equation  of  the 
form 


m  6  +  kS  *  -P$t  +  f  on  1*^  (4.2) 

where  P  is  the  (assumed  constant  in  linear  theory)  density  of  the  fluid, 
m,k  are  positive  constants,  and  f  is  an  applied  external  force  on  r  .  if 
we  assume  the  boundary  section  r  is  impenetrable  we  obtain  a  third  equation 
from  the  continuity  of  velocity  at  the  boundary, 

't 

3$ 


5  =  -r1  on  T 

t  on  1 


(4.3) 


where 


3n 


is  the  outward  normal  velocity. 


(This  model  has  been  given  in  [5],  p.  263;  a  mathematical  analysis  is  found  in 


(This  model  has  been  given  in  [5],  p.  263;  a  mathematical  analysis  is  found  in 


[6],  and  a  discussion  of  linear  control  synthesis  appears  in  [7].) 

If  we  assume  G  is  bounded  and  F^  u  u  *  3g  where  is  a  free 

surface  then  the  pressure  on  must  take  on  the  ambient  atmospheric  pres¬ 

sure  distribution.  If  we  normalize  this  boundary  pressure  distribution  to 
zero  we  find 

=  0  on  r2>  (4.4) 

r3  is  a  rigid  boundary  so 

-  0  on  T3.  (4.5) 

In  this  paper  we  consider  the  case  where  G  is  a  right  circular  cylinder 
with  free  surface  at  the  top  and  reacting  boundary  at  the  bottom.  This  is 
shown  in  Figure  2. 


3£ 

3n 


x 


Figure  2 
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For  simplicity  w«  set  c-P-m-k-  1.  Also  we  seek  a  solution  of 


(4.1)  -  (4.5)  ♦  •  f(x,t).  In  this  case  (4.5)  is  automatically  satisfied. 

Differentiation  of  (4.1)  -  (4.4)  with  respect  to  t  yields  the  system 


wtt  "  wxx  0  <  x  <  1 

®tt  +  ®  *  -wt(t,1)  +  v(t)w(t,1)#  X  -  1, 

9  “  w  ,  x  «  1 , 

t  X 

w  *  0  ,  x  ■  0, 


(4.1') 

(4.2*) 

(4.3*) 

(4.4*) 


where  8  5  6  and  we  have  synthesized  the  external  driving  force  f ' (t)  = 
v(t)w(t,1),  i.e.,  we  are  looking  for  a  feedback  which  then  multiplies  the 
observed  pressure  at  x  *  1 .  of  course  we  must  also  specify  initial 
conditions 

w(0,x)  -  w0(x),  wt(0,x)  -  w1 (x) ,  9(0)  -  8q,  0^0)  m  91 '  0  <  x  <  1. 
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5.  Abstract  fromulation  of  hybrid  system  2  and  feedback  stabilization. 
For  the  hybrid  system  (4.1')  -  (4. 4')  set 


r 

w 

'  r. 

w  > 

t 

0 

w 

w 

0 

t 

8 

,  Au  = 

XX 

8 

,  8u  = 

0 

t 

6tJ 

-8  -  wfc ( • , 1 )  j 

w{  •  ,  1  )  j 

X  =  {(w,w  ,0,6  ).*  w  e  ^(0,1),  wt  e  L2(0,1),  8  e  R,  0t  e  R,  w  =  0  at  x  =  0} 
endowed  with  inner  product 

*  ,  i  *  *  *  * 

<u,u  >  =  J  (w  w  +  w.w  )dx  +88  +  98  , 

XqXX  tt  tt 

d(  A)  =  {(w,tt,8,e  )  e  x;  we  h2(o,d,  wfc  e  h^o.d, 

w  =0  at  x  =  0  ,  8  =w  at  x=1 }. 
t  t  x 

Since  for  w  e  H^(0,1)  the  Sobolev  lemma  ([4],  p.  174])  says  w  can  be  viewed 

as  an  element  of  C(0,1)  after  possible  modification  on  a  set  measure  zero, 

we  shall  take  such  w  as  a  continuous.  So  for  w  e  H'(0,1),  w  =  0,  w  is 

identified  as  a  continuous  function  which  satisfies 

w(x)  =  /x  w' (x}dx 
0 

in  the  Lebesgue  sense.  Hence  the  boundary  value  w(l)  is  simply 

w(1)  =  /1w'(x)dx.  (5.1) 

0 

It  is  in  this  sense  that  the  boundary  conditions  in  X  and  D(A)  are  to  be 
understood. 

We  note  for  u  e  D(A),  <Au,u>v  =0  so  A  is  certainly  dissipative.  A 
also  densely  defined  and  Range  (I-A)  =  X  (by  a  direct  computation)  so  the 
Lumer-Phillips  Theorem  [4]  implies,  A  is  the  infinitesimal  generator  of  a 
Cq  semigroup  of  contractions  e^fc  on  X.  Also  we  note  that  if  {u^}  is  a 
bounded  sequence  in  X  ,  (5.1)  shows  {w^ ( 1 ) }  to  be  a  bounded  sequence  in 
R.  Hence  it  possesses  a  convergent  subsequence.  Thus  8  :  X  +  X  is  compact. 


is 
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All  the  hypotheses  of  Theorem  1  are  satisfied  modulo  showing  that  the  only 
At  At  + 

solution  of  <e  e  i(>>  =*  0  for  all  t  e  R  is  =  0.  In  this  problem 

<eAtV,  eAt<l»  =  9tw(1,t) 

where  w,8  are  solutions  to  the  uncontrolled  (p  =  0)  system  (4.1* >  - 

(4.41).  Separation  of  variables  shows 

~  *  io  t  _  -io  t 

w(x,t)  =*  l  (A^e  +  A^e  K  }  sin  a^x  ,  (5.2) 


where  a 

Here 

equation 


_  -  ifl  t  _  -io  t 

®t(t)  -  l  (Ake  +  A^e  )a^cos<\ 
k“1 

are  the  positive  roots  of 


k  k 


tana  - - a. 

a 

denotes  complex  conjugate.  If  we  substitute  (5.2), 
Q^wd/t)  **  0  for  all  t  e  R+,  we  obtain 

2  l  |Ak|2aksin  a{  ~os  ay 
k-1 


•  _  i20^t  -i20  t 

+  l  (\e  +  \  ®  *  )°k  sin  \  COB  °k 

k-1 


■  i(o  +a  )t  i(a  -a  )t 

,  v  .  nk  t  -  -  nk 

+  1  (A  A  e  +  A  A  e 

,  .  k  n  n  k 

n,k*1 

n*k 


(5.3) 


(5.4) 

(5.3)  into  the 


i(°n"°k)t  -  -  i(ak"an,t) 

+  A  A.  e  +  A  A.  e  a  sin  a  cos  a 

n  k  n  k  n  k  n 

-  0  for  all  t  e  R+.  (5.5) 

The  terms  on  the  left  hand  side  of  (5.5)  form  an  almost  periodic  function  in 

t.  If  none  of  the  frequencies  appearing  in  the  second  sum  appear  in  the  first 

and  third  sums  the  uniqueness  theorem  for  almost  periodic  functions  will  imply 

Ak  “  0,  k  *  1,2, •••  .  This  in  turn  will  imply  w(x,t)  *  0,  9  (t)  =  0  which 

*  t 

At  At 

with  (4.2')  shows  =  0  will  be  the  only  solution  of  <e  w,  e  <J>>  =  0  for 

all  t  e  R* .  So  we  now  prove  the  following  lemma. 
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Lemma  7.  Let  o  be  the  positive  roots  of  (5.4),  0  <  o„  <  o^  <  o„  • • • 

"  n  \  Z  3 

Then  a  ±  a  ~  20,  holds  only  when  r=s=k  and  +  sign  is  taken, 
r  s  k 

Proof.  For  later  use  we  record  the  first  three  positive  roots  of  (5.4): 

0i  =  0.67625,,» 

02  =  2. 1 1708* •• 

03  =  4.92125***  . 


Case  1 :  Assume  o^  +  =  2 0^.  Then  taking  tan  of  both  sides  and  using  (5.4) 


we  find 


1-00 
_ r  s 


40^-30  o  -(1-o  O  )' 
k  r  s  r  s 


Solving  for  o  o  we  see 
r  s 


o  o 
r  s 


(5.6) 


2  2 
(Note  o  o  *1/0*1  so  our  manipulations  are  valid).  If  o  <  2/3  or 
r  s  ic  k 

2 

Ok  >  1  we  see  arag  <  0  which  violates  positivity  of  the  roots.  Since 
2 

o  (2/3,1)  for  all  k,  we  see  that  0  +  0  =  2 0,  holds  only  when 

k  r  s  k 

o  =  o  -  0  ,  i.e.,  r=*s=k. 

r  s  k 


Case  2.  Assume  ar  ”  °a  =  2a^  and  proceed  as  in  Case  1.  We  find  this  time 
that 


„2  ,  ,  1 
o  =»  1  +  -  , 

k  oo-3 

r  s 


(5.7) 


If  r,8  >2,  r  *  s,  graphical  analysis  of  (5.4)  shows  o  -  o  >  o  -  o 

r  s  3  2 

2 

*  2.80417***  .  So  if  20^  “  ar  ~  ag  we  must  have  0^  >  1.965#*»  .  Also 

we  know  for  r,s  >2  o  o  >3.  Hence  from  (5.7)  we  see  that  if 

r  s 


0  -  0  -  20  for  r,s  >  2,  r  *  s,  that 

IT  3  JC 
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r 


1.965* • •  <  1  +  - 

0  0-3 
r  s 

for  3  <  °r°s  <  4.035* ••  .  No  pair  ar>ag  satisfy  this  relation. 

On  the  other  hand  for  r  or  s  <  2  we  can  only  have  s  =  1 ,  r>1.  We 

consider  the  set  {o  }  so  that 

r 

1  +  g~o'-j"  <  °2  “  4.48203***  .  (5.8) 

1  r 

This  set  is  the  set  of  {o^}  80  that  °iar  >  3.287186***  or  o ^  >  4.8609***  , 

i.  e.,  r  >  3.  So  for  r  >  3  the  only  way  (5.7)  can  hold  with  s  =*  1  is 

when  k  ■  1.  i.e.,  o  -  o  «  20„  or  o  *  30  .  But  there  is  no  such  o  . 

r  1  1  r  1  r 

Finally  if  r  ■  2  we  must  have  0  -  0^  ■  2°k'  *,e*'  ak  *  *720415  which 

also  cannot  hold.  The  lemma  is  proven.  U 

We  can  now  conclude  from  the  discussion  proceeding  Lemma  7  that  the 
following  theorem  is  true. 

Theorem  4.  The  feedback  system  (4.1*)  -  (4.4*)  with 

v(t)  -  -  et<t)w(i,t) 

possesses  a  unique  globally  defined  weak  solution  for  initial  data 
(w  ,w  ,8  ,8  )  e  X  and  (w  ,  w  ,0  ,0  )  ♦  (0,0, 0,0)  weakly  in  X  as  t  ♦  • 

vJ  1  U  1  t  v 

and  (4.1  * )-(4.4’ )  is  weakly  stabilizable. 
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